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Abstract 



We investigate in stable strange hadronic matter (SHM) the modifica- 
tion of the masses of the scalar (a, a*) and the vector (u, (ft) mesons. The 
baryon ground state is treated in the relativistic Hartree approximation in 
the nonlinear a-uj and linear a*-(ft model. In stable SHM, the masses of all 
the mesons reveal considerable reduction due to large vacuum polarization 
contribution from the hyperons and small density dependent effects caused 
by larger binding. 

PACS: 21.65+f, 24.10Jv 
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The study of the properties of hadrons, in particular the light vector mesons (u, p, 0), has 
recently attracted wide interest both experimentally and theoretically. Recent experiments 
from the HELIOS-3 |1| and the CERES JJ] collaborations indicate a significant amount of 
strength below the p-meson peak. This has been interpreted [||] as the decrease of the p- 
meson mass in the medium. In an effective chiral model based on the symmetries of QCD, 
Brown and Rho [[| predicted an approximate scaling law for the in-medium decrease of the 
masses of the mesons and nucleons. On the other hand, in the quantum hadrodynamic 
(QHD) model |5],|6|] based on structureless baryons, the meson masses are found to increase 
when particle-hole excitations from the nucleon Fermi sea are considered. By also including 
the particle-antiparticle excitations from the Dirac sea, the meson masses has been, however, 
shown to decrease |7|-f§. 

So far all investigations of in-medium meson mass modification has been done in the 
nuclear matter environment. However, presently there is a growing interest about the possi- 
bility of bound strange matter. In analogy to the stable strange quark matter [JTUJ , strange 
hadronic matter (SHM) composed of nucleons and hyperons has been speculated to be 
absolutely stable or at least metastable with respect to weak hadronic decays. It is expected 
that such strange matter may possibly be created at RHIC and LHC [0. It is therefore 
worth investigating the properties of meson masses in the strange baryonic matter environ- 
ment. This study is not only of interest in itself but could serve as a signal of the formation 
of SHM in heavy ion collisions. In this letter, in a complete self-consistent calculation within 
the framework of the QHD model, we examine the mass modifications of both, the scalar 
and the vector mesons in SHM. 

Let us consider the composition of SHM which is stable with respect to particle emission. 
The analysis of level shifts and widths of S~ atomic level suggest [H| a well-depth of E 
in nuclear matter of ~ 20 — 30 MeV. While for strong processes, EiV — > AN and 

SA — > EN, the energy released is Q ~ 78 and 52 MeV, respectively. Consequently, systems 
involving E's are unstable with respect to these strong decays. Analysis of binding energies of 
A, and emulsion experiments with K~ beams [TB[ yields well-depths of A and H in nuclear 



matter of U^Q ~ 28 — 30 MeV. However, in strong decays, AA ^ NS, the system can 
be stable since Q ~ 25 MeV. Therefore, we consider here only the set of baryon species 
B = {N, A,H}, which can constitute stable SHM. 

Up to now all investigations of meson mass modification has been performed in the simple 
linear Walecka model with nucleons only. Considering here nonlinear self-interactions of the 
scalar field a in the QHD model, the total Lagrangian is 

£ = (*7^/i -m B + 9aBO- - 9uBln^) lp B 

B 

+ l - (d,ad»a - mla 2 ) + U(a) - \uj^ v + \mlu^ + 6C , (1) 

where the summation is over all baryon species B = {N, A,S}. The scalar self-interaction 
U(a) = 5 , 2 " 3 /3 + (73<r 4 /4 yields sufficient flexibility to the effective Lagrangian of the model 
and they are also necessary for a good reproduction of nuclear ground state properties, in 
particular the compressibility. The term 8L contains renormalization counterterms. This 
model is not able to reproduce the observed strongly attractive AA interaction. The situation 
can be remedied by introducing two additional meson fields, the scalar meson / (975) (de- 
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noted as a* hereafter) and the vector meson 0(1020) [11] which couple only to the hyperons 
(Y). The corresponding (linear) Lagrangian is 
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The self-consistent propagator in the medium for baryon species B can be written as sum 
of the Feynman [G^(A;)] and density-dependent [G^(fc)] parts: 



G B {k) = Gl{k)+G%(k) 



+ 
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k* 2 -m*^ + ir] ' E* B {k) 



6(k* -E B (k))9(k FB -\k\) 
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The momentum and energy of baryon B are k* = (ko — g^B^o — <?<£B0o,k) and E* B {k) = 
[k 2 + m^ 2 ] 1//2 , and kp B are the Fermi momenta. The scalar mean fields shift the mass mg of 
the baryons both in the Fermi and Dirac sea to 
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In mean field theory (MFT), the total energy of the system is generated by the presence of 
all baryons in the occupied Fermi seas. In contrast, in the relativistic Hartree approximation 
(RHA), the effect of the infinite Dirac sea is also included. The renormalized total energy 
density in RHA is given by 
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where £mft is the usual MFT energy density [TTJ . The contribution from vacuum fluctuation 
of all the baryons is given by 
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where I B = 21+ 1 is the isospin degeneracy of the baryon B. The mass shift of the Dirac sea 
from m B to m B produces the vacuum fluctuation contribution. The renormalized nonlinear 
cr-meson contribution is [HI 
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1 + A x + A 2 ) 2 ln(l + Ai + A 2 ) - (Ai + A 2 ) 



-^(Ai + A 2 ) 2 -^A 2 (A 1 + 3A 2 ) + lAf 



(7) 



where Ai = 2giojm\ and A 2 = Sg^a 2 /m 2 a . 

The meson propagator in the baryonic medium can be computed by summing over bub- 
bles which consists of repeated insertions of the lowest order one-loop proper polarization 
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part. This is equivalent to relativistic random phase approximation (RPA). Since both, 
scalar (a, a*) and vector (u>, 0) mesons are present in our model, it is essential to include 
scalar- vector mixing which is a pure density dependent effect. Therefore it is convenient 
to define a full scalar-vector meson propagator in the form of a 5 x 5 matrix with 
indices a, h ranging from to 4, where 4 corresponds to the scalar meson and to 3 the 
components of vector meson. Moreover, the (strangeness violating) scalar-vector coupling 
between the strange and nonstrange mesons are prohibited by invoking the OZI rule. We 
therefore consider couplings between <j-uj and a*-<f) separately. We shall present explicitly 
the calculations only for the a-u propagator; the expressions for a*-(f> propagator would 
follow similarly. Dyson's equation for the full a-uj propagator T> can be written in matrix 
form as 

V = V° + V°UV , (8) 
where T>° is the lowest order a-u meson propagator 

v° = ( D r A °oV o) 



A ( . 

expressed in terms of the noninteracting a and tu-meson propagators 

A°(ff) = 2 * . , (10) 

% ~ m a + IV 

D%(q) = (g, u -q,q„/ml)D°(q) , (11) 

D°(q) = 2 ~* . , (12) 
% -mi + 17] 

where q 2 = q$ — q is the four-momentum carried by the meson. Note that the effect of 
nonlinear a-interaction (the boson loops) is to replace the bare meson mass m 2 a in Eq. (10) 



by ml = ml + d 2 U (a) / da 2 = m 2 a + 2g 2 a + 3g 3 a 2 (see Ref. JT4 



The polarization insertion of Eq. (8) is also a 5 x 5 matrix 

Ub ni M)B (9) K(i) + Eb (g) ) ' [L6) 

where each entries in Eq. (13) is summed over all the baryons, except for the pion loop H*(q). 
The effect of relatively large cr-width has been accounted by including the contribution of 
pion loop to U a [§]; the in-medium modification of the pion loop is neglected. The pion 
propagator A v is obtained from Eq. (10) with m 2 replaced by m 2 . The pion loop polarization 
contribution to the a-meson is 

K(<l) = \w™™l J -0^ A n (k)A n (k + q) . (14) 

In terms of the baryon propagator the lowest order a, u and a-u (mixed) polarizations for 
the baryon loop B are respectively given by 
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nf )B (g) = ig„ B 9.B / 7^1 Tr [G B (k) lfl G B (k + q)} . 
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As in the case for the baryon propagator (Eq. (3)) the above polarization insertions (except 
ITJ) can be expressed as the sum of Feynman (F) part and density-dependent (D) part, i.e. 

U B = U B(F) + U B(D)_ finite D _ part hag form qD x qD _|_ qD x qF _|_ qF x qD 

describes particle-hole excitations and also includes the Pauli blocking of BB excitations. 



Since the polarizations of the D-part are defined in Ref. ||15|| , we will not refer them here. 

The divergent F-part of the polarization insertions can be rendered finite by adding 
appropriate counterterms to the Lagrangian of Eq. (1). For the a, each of the baryons loops 
and the 7r-loop have to be renormalized separately. For any baryon loop contribution to the 
a the usual counterterm Lagrangian |5||J is used 
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The coefficients a<i and ( a can be obtained by imposing the condition that the propagator 
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in vacuum (m* B = m B ) reproduces the "physical" properties of a-meson |~5 
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The renormalized cr-meson self-energy for a baryon loop is 
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For the 7r-loop we employ the renormalization condition in free space, Il^^ RF \q 2 ) = at 



m l 0- We obtain finally 

K {RF) (a) 
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For the lu, only a wavefunction counterterm = ^uj^uj^ /A is required to make the 
F-part Il£, finite. The renormalized ui self-energy is Ti^ l } RF " > = {—g^q 2 + q^q u )H B ^ RF \ 
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Employing the renormalization condition in vacuum, Tl B ^ RF \q 2 ; m* B = tub) = at q 2 = m^, 
we obtain 

As mentioned before, the mixed part of Eq. (17) does not contribute to vacuum 

polarization. 

The solution of Dyson's equation, Eq. (8), is T> = P°/(l — D°TT). By defining the 
dielectric function e as 

e = det (l - X?°n) = 4 e L , (23) 

the poles of the a-u propagator which define the respective meson masses in the medium are 
now contained in the dielectric function when e — 0. By taking q = (0, 0, q) where q — |q|, 
we obtained in Eq. (23) the transverse and longitudinal dielectric functions defined as 

e T =(l + D°U T ) , (24) 



e L =(l- A°n CT ) (l - D°IL L ) + | A°D° (nJ M) ) 2 , (25) 

where the polarization insertion of Eq. (16) is now split into transverse Ut = Tin = TI22 
and longitudinal TX L = Tl 00 — IT33 components. Because of baryon current conservation, 
only the 0-th component of the mixed part Tl M survives. Note that in Eqs. (24) and (25) 
the polarizations represent those which have been summed over all the baryons (see Eq. 
(13)). The eigencondition for determining the collective excitation spectrum (i.e. finding 
the effective meson masses) is equivalent to searching for the zeros of the dielectric function. 
In particular, for a given three-momentum transfer q = |q|, the "invariant mass" of a meson 

(a or uj) is m* m = \Jq$ — q 2 , where q is obtained from the condition e — 0. In the present 
study of meson mass modifications in the medium, we restrict ourselves to the meson branch 
in the time-like region (q 2 > 0). 

Since the propagation of the strange a*-(ft mesons are decoupled from that of the non- 
strange cr-u mesons, we may follow the same procedure as given in Eqs. (8)- (25) to obtain 
the effective masses of the strange mesons. In particular, a dielectric function (Eq. (23)) is 
obtained but with the masses of the mesons and their couplings to the baryons correspond 
to the a* and <fi. Moreover, since a linear a* interaction is used in this case, m 2 a in Eq. (10) 
should be replaced by the bare meson mass m 2 a » . 

The renormalization conditions in RHA in Eq. (6) is imposed at q 2 = 0, while those used 
in Eq. (19) for a and cr*-mesons are at q 2 = m 2 ^ and q 2 = ml*, respectively. This difference 
yields an additional term in the vacuum fluctuation energy [T5j , so the total energy density 
for SHM is 

£ = £rha + ^ ^ T B !b m B( m *B ~ m B? , (26) 

B 47r 

where £rha is the RHA energy of Eq. (5) and 
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a aB = -. 9- + - / 

Am B 2 Jo 

and a similar expression for a CT *£ for a*; for nucleons a CT »£ = 0. 

The field equations are obtained by minimizing the energy density £ with respect to that 
field. At a given baryon density % and strangeness fraction fs = (n\ + 2ri=)/ns, the set 
of field equations are solved self-consistently in conjunction with the chemical equilibrium 
condition 2fi\ = fi^ + due to the reaction AA ^ NE. The chemical potential of a baryon 
species B is hb = [^f s + rn* B 2 ] 1 / 2 + g^B^o + 9<f>B4>o- The four saturation properties of nuclear 
matter (NM): density n = 0.16 fm -3 , binding energy E/B = —16 MeV, effective nucleon 
mass m* N /rriN = 0.78 and compression modulus K = 300 MeV are used to fix the nucleon 
coupling constants g a n, Q^n and the parameters #2 and #3 of the a self- interaction. The 
coupling constants for pure NM without hyperons are shown in Table I. When hyperons are 
included, i.e. for SHM, they will contribute to S from their vacuum fluctuations Vb even 
if their Fermi states are empty. This entails a redetermination of the coupling constants 
for the nucleons. The Vb depends on the effective baryon mass, which in turn depends 
on the scalar-baryon coupling constants (see Eqs. (4) and (6)). Therefore, the a and a* 
couplings to the hyperons (Y) should be predetermined. For this purpose, we adopt the 
SU(6) model, i.e. g a A/g aN = 2/3, g aS /gaN = 1/3 for the a-Y couplings, and g a *klg<,N = 
a/2/3, ga's/gaN = 2a/2/3 for the a*-Y couplings. Note that the nucleons do not couple to 
strange mesons, i.e. g a *N = g^N = 0. The coupling constants of nucleons for the SHM are 
given in Table I. The couplings of u to the hyperons can be obtained from the well-depth 
of Y in saturated NM: Uf ] = g aY v eq - g^ eq with U ( A N) = 30 MeV and U^ N) = 28 MeV. 
The — Y couplings are obtained by fitting them to a well-depth C/^a ~ 40 MeV, for a 
A or H in a H "bath" with n= ~ n [|TTJ. To determine the couplings of the pion to scalar 
mesons, g a7r and g a *ir, we adjust them to reproduce the widths of a and a* in free space, i.e. 
T° = —^sU^ RF ^ /m s at q 2 = m 2 s . For a conservative estimate we consider T° = 300 MeV 
and = 70 MeV. 

The stability of strange hadronic matter may be explored by considering its binding 
energy defined as E/B = S/ns — J2i Y%m>i, where the abundance Yj, = rii/nB- In Fig. 1, we 
present the binding energy E/B as a function of baryon density at various strangeness 
fractions f$- With increasing f$, the binding energy of SHM is found to increase and the 
saturation point is shifted to higher density. This is a consequence of the opening of new 
(strangeness) degrees of freedom in the form of hyperons. At high densities for large fs 

TABLE I. The nucleon-meson coupling constants in RHA obtained by reproducing the nuclear 
matter saturation properties (see text). The results are for nuclear matter (NM) with strangeness 
fraction fs = and for strange hadronic matter (SHM) which includes also the vacuum polarization 
of hyperons. All the couplings are dimensionless except 52 which has the dimension of fm" 1 . The 
bare hadronic masses are = 939 MeV, m\ = 1116 MeV, m= = 1313 MeV, m a = 550 MeV, 
m u = 783 MeV, m CT * = 975 MeV, and = 1020 MeV. 
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values when the Fermi energy of nucleons exceeds the effective masses of the hyperons minus 
their associated interaction, the conversion of nucleons at its Fermi surface to hyperons of 
increasingly larger mass is favored. This conversion lowers the Fermi energy of the nucleons 
leading to enhanced binding. A stable SHM with a maximum binding of E/B = —38.46 
MeV is obtained for a large fs ~ 1.35 at a relatively high density ub ~ 3uq. A further 
increase of fs enforces the Fermi energy of the hyperons to increase resulting in the decrease 
of binding. This finding is quite similar to that obtained in the quark-meson coupling model 

The variation of baryon effective masses m* B is shown in Fig. 2 as a function of density 
jig for varying strangeness fs- It is observed that m* B decreases with increasing n B , with 
m* N having the largest decrease rate and m| the smallest at each fs value. Furthermore, 
with increasing fs the effective nucleon mass increases while the effective masses of hyperons 
decrease at any density. This effect stems from decrease of the nonstrange meson fields a and 
lo and increase of the strange meson fields a* and with increasing fs- Since the nucleons 
couple only to a and u>, m* N is increased. The hyperons, however, do couple to all the meson 
fields resulting in a decrease of their effective masses with increasing fs, especially for H 
which has the strongest coupling to the strange mesons. 

In Fig. 3, we display the "invariant mass" m* m = \J q$ — q 2 of nonstrange a-meson, 
m*, (left panels) and cu-meson, m* , (right panels) as a function of baryon density % for 
different strangeness fractions fs- We consider first the features observed for m* for small 
three-momentum transfer q = |q| = 1 MeV (top- left panel). For fs = 0, m* is found to 
decrease with increasing % for small values of n B < n . This reduction is caused due to two 
competing effects. The vacuum polarization which leads to reduction of m* dominates over 
the density dependent dressing of the meson propagator which causes an increase in m*. In 
fact, this reduction can be traced back to the corresponding reduction of m* N < uin in the 
medium (see Eq. (20)). However, the decrease in the scalar polarization depends on m* B , 
m B , and m 2 in a complicated fashion. For densities above n , the density-dependent (D) 
part becomes increasingly dominant resulting in increase of m*. In SHM with fs ^ 0, the 
vacuum polarization contribution from the hyperons and the nucleons causes a considerable 
suppression of m* at low densities. At higher densities m* increases, however, the rate of 
increase above ~ no becomes smaller with increasing fs- An explanation to this effect is 
as follows. The D-part of a propagator is primarily determined by nucleons to which it 
has the strongest coupling. The increase of fs causes an increase in the Fermi momenta 
of the hyperons while that of the nucleons decrease. Consequently, with increasing fs the 
decreasing D-part of the nucleons results in a slower rate of increase of m* with ris- 

At small value of q — 1 MeV, the (density dependent) effect of scalar-vector mixing is 
negligible (see Eq. (25)). On the other hand, for large values of q = 500 MeV (central-left 
panel) and q = 1 GeV (bottom-left panel), the mixing is more effective. (For the latter 
value, the present model may have been stretched to the extremes of applicability.) It gives 
rise to a repulsion which again leads to decrease of m* at high density. With increasing fs, 
the onset of this decrease or the peak positions in the figure are found to shift to higher 
densities. This is a manifestation of the shift of saturation value E/B to higher ub as fs 
increases (see Fig. 1). 

For small values of q = 1 MeV, the transverse and longitudinal invariant cu-meson mass, 
m* T and m* WL are practically identical as evident in Fig. 3 (top-right panel). In contrast 
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to cr mass, the vacuum polarization effect is much stronger (see Eq. (22)) and the density 
dependent effect is much weaker for the heavier u mass. This causes for f$ = a substantial 
decrease of m* up to n B rs 2n , and a subsequent small increase with density. The reduction 
is more pronounced for SHM where the total vacuum polarization effect is stronger and the 
D-part (mainly controlled by N) is weaker. In fact, for f$ = 1.5, decreases considerably 
even up to large densities. As a consequence, the crossing between m* and m* at q — 1 
MeV is shifted to lower ub for the SHM. 

At large values of q = 500 MeV and 1 GeV, the longitudinal mass m* L (solid line) and 
transverse mass m* L (dashed line) get well separated. It is found that m* L is reduced near 
nuclear matter density and finally increases with density attaining values higher than m* . 
As in the q = 1 MeV case, the reduction in mass is stronger for SHM. Because of strong 
repulsion from the mixing at these high three-momentum values, m* and m*^ never cross 
each other. 

In Fig. 4 the "invariant mass" of the strange <r*-meson, m* » , (left panels) and 0-meson, 
m*i, (right panels) are shown as a function of baryon density % for different strangeness 
fractions. In this case the masses are determined by the hyperons as nucleons do not couple 
to the strange mesons. As for the nonstrange mesons, the masses m*» and rri^ in general 
decreases with increasing ns- However, the decrease is more enhanced over the large density 
range explored here. This arises due to large vacuum fluctuation contribution from the 
hyperons, and, in particular, a small density- dependent part of the hyperons stemming from 
the large binding in SHM. At high densities % > 4no, in contrast to nonstrange meson 
masses, the strange meson masses has higher values for larger f$. This reversal in behavior 
for the m** and results from the increase of the Fermi momenta and hence the D-part 
of the hyperons for large f$- The scalar mass m** however becomes sensitive to high values 
of q and drops at fs = 1.5 below that for fs = 0.5. 

We have investigated the meson mass modification in strange hadronic matter as a 
function of baryon density for a fixed strangeness fraction. The results correspond to a 
metastable SHM. For a given ns, an absolute stable SHM could be obtained by determining 
the absolute minimum of the binding energy E/B as a function of fs- In this situation, 
we have found that the masses of the mesons in SHM undergo a drastic reduction over 
a wide range of density. Besides the large vacuum part, the minimum in the energy per 
baryon E/B at each density enforces a smallest possible Fermi momentum and hence the 
density-dependent for all the baryons. 

It is worth mentioning here that in the linear Walecka model, it has been demonstrated 
PJ6|J17[| that a self-consistent inclusion of exchange (Fock) term has an insignificant contribu- 
tion to binding energy when the two parameters of this model, g a ^ and g^N, are renormal- 
ized to reproduce the same NM saturation properties of baryon density and binding energy. 
Even the predicted values of the effective nucleon mass, m* N , and incompressibility, K, in 
this Hartree-Fock calculation are almost identical to the Hartree results ||. However, the 
applicability of the linear Walecka model to moderate and high density phenomena can be 
misleading if the two parameters m* N and K at n Q are not under control. In fact, it was 
shown [[DJ that in the nonlinear Walecka model when the four parameters, g a N, g^N, <?2, 
and #3 (the latter two are from the a self-interaction) are fitted to the same NM saturation 
properties of n , E/B, m* N , and K (as in our calculation), nearly all the properties obtained 
in the relativistic Hartree calculation differ from the mean field results only by m 3% even at 
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n = IOuq. (This is in contrast to the linear Walecka model results.) It is therefore expected 
that by further inclusion of exchange terms (in all the baryonic sectors) and pseudoscalar 
mesons in the present nonlinear Walecka model, and performing a self-consistent calcula- 
tion (with the parameters renormalized to the same NM saturation properties), the results 
for E/B and m* B will be practically unaltered from the Hartree calculation. Consequently, 
exchange corrections from the nucleonic and strangeness sectors should also have an insignif- 
icant contribution and therefore neglected in the present study of the modification of meson 
masses in the medium. 

In summary, we have investigated the masses of baryons (N, A, S) and, in particular, 
the nonstrange (a, uj) and strange (a*, 0) mesons in stable strange hadronic matter. The 
ground state properties of the SHM in the relativistic Hartree approximation is obtained by 
using a nonlinear a-uo and linear a*-<p Lagrangians. With increasing strangeness fraction, 
the effective mass of the nucleons increases while that of the hyperons decreases. The masses 
of all the mesons reveal a considerable reduction over a wide density range with increasing 
strangeness. This may be attributed to a large contribution from the vacuum polarization of 
the hyperons which causes the decrease in the meson masses at small densities. The larger 
binding for the SHM and therefore a smaller density-dependent part helps to reduce meson 
masses at high densities. 



S.P. and S.G. acknowledge support from the Alexander von Humboldt Foundation. 
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FIG. 3. The in-medium masses of a-meson (left panels) and u>meson (right panels) as a function 
of baryon density for various fs with q = |q| = 1, 500, and 1000 MeV. The a- meson masses are 



14 



for fs from to 1.5 in steps of 0.5. In the right panels, the solid and dashed lines correspond to 
longitudinal and transverse w-meson masses, respectively. 
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FIG. 4. Same as Fig. 3 but for the cr*-meson (left panels) and c/>-meson (right panels). 
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